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Abstract 

Two fundamental objects in knot theory are the minimal 
genus surface and the least area surface bounded by a 
knot in a 3 -dimensional manifold. When the knot is em- 
bedded in a general 3-manifold, the problems of finding 
these surfaces were shown to be NP-complete and NP- 
hard respectively. However, there is evidence that the spe- 
cial case when the ambient manifold is U 3 , or more gener- 
ally when the second homology is trivial, should be con- 
siderably more tractable. Indeed, we show here that a nat- 
ural discrete version of the least area surface can be found 
in polynomial time. 

The precise setting is that the knot is a 1 -dimensional 
subcomplex of a triangulation of the ambient 3-manifold. 
The main tool we use is a linear programming formulation 
of the Optimal Bounding Chain Problem (OBCP), where 
one is required to find the smallest norm chain with a 
given boundary. While the decision variant of OBCP is NP- 
complete in general, we give conditions under which it 
can be solved in polynomial time. We then show that the 
least area surface can be constructed from the optimal 
bounding chain using a standard desingularization argu- 
ment from 3-dimensional topology. 

We also prove that the related Optimal Homologous 
Chain Problem is NP-complete for homology with integer 
coefficients, complementing the corresponding result of 
Chen and Freedman for mod 2 homology. 

1. Introduction 

A knot K is a simple closed loop in an ambient 3-dimen- 
sional manifold Y. Provided K is null-homologous, which 
is always the case if Y — R 3 , there is an embedded ori- 
entable surface S in Y whose boundary is K (equivalently 
S is a compact smooth orientable surface in Y without 
self-intersections and with boundary K). A fundamental 
property of K is the minimal genus of such an S, which is 
denoted g[K] (we take g{K) = oo if there are no such sur- 
faces). In the 1960s, Haken used normal surface theory to 
give an algorithm for computing g{K], opening the door 
to a whole subfield of low- dimensional topology and lead- 
ing to the discovery of algorithms for determining a wide 
range of topological properties of 3-manifolds [18]. How- 



ever, algorithms based on normal surface theory are quite 
slow in practice [2-4], and there are very few results that 
have been verified via such normal surface computations 
[5, 10]. Moreover, in some cases the underlying problems 
have been shown to be fundamentally difficult. In their 
foundational work, Agol, Hass, and Thurston showed that 
the following decision problem is NP-complete [1]: 

1 . 1 Knot Genus. Given an integer go and a knot K embed- 
ded in the 1-skeleton of a triangulation of a closed 3-mani- 
fold Y, is g{K) < go? 

While Knot Genus is NP-complete, when Y is orientable 
and the second Betti number b 2 (Y) = rank(_ff2(F; Z)) is 0, 
for instance Y = S 3 , then this problem likely simplifies. 
While their project is not yet complete, Agol, Hass, and 
Thurston have developed a very promising approach to 
showing that when b 2 {Y) = there is a certificate for the 
complementary problem g(X) > go which can be verified 
in polynomial time. This would mean that this special case 
of Knot Genus is also in co-NP, raising the possibility of 
a polynomial-time algorithm when b 2 [Y) - 0. However, 
currently there are no known algorithms which exploit the 
fact that b 2 {Y)- 0. Despite this, our long-term goal is 

1.2 Conjecture. For orientable Y with b 2 {Y) = the Knot 
Genus problem is in P. 

Here, we study the related problem of finding the least 
area surface bounded by a knot. This problem has its ori- 
gin in classical differential geometry, as we now sketch 
starting with the case where the ambient manifold Y is R 3 . 
For a smooth knot K in U 3 there is always a smooth embed- 
ded orientable surface ScR 3 with dS = K. By deep the- 
orems in Geometric Measure Theory, there always exists 
such a surface So of least area [19]. A least area surface So 
is necessarily minimal in that it has mean curvature ev- 
erywhere, like the surface of a soap-bubble. It is typically 
impossible to find the least area surface analytically, and 
the first paper on numerical methods for approximating 
So appeared in 1927 [9] . An algorithm to deal with arbitrary 
K was first given by Sullivan [27] in 1990; see also [21-25] 
for alternate approaches and numerical experiments. 

Of course, one can consider this question for null- 
homologous knots in an arbitrary Riemannian 3-mani- 
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fold Y, and one has the same existence theorems for least 
area surfaces when Y is closed. Agol, Hass, and Thurston 
considered a certain discrete version of this problem, and 
showed that the question of whether K bounds a surface 
of area < Aq is NP-hard [1] . Because they put no restriction 
on the surface involved, it is not clear if their question is 
in NP. Here, we consider another discretization, the Least 
Spanning Area Problem of Section 5, which is a little more 
combinatorial and will thus turn out to be NP-complete 
(Theorem 5.2). For this problem, we show 

1.3 Theorem. For orientable Y with bziY) = the Least 
Spanning Area Problem is in P. 

In Section 6, we discuss how the ideas behind Theorem 1.3 
might be used to attack Conjecture 1.2, as these two ques- 
tions have a very similar flavor. 

One of two key tools behind Theorem 1.3 is the fol- 
lowing type of combinatorial optimization problem. For 
a finite simplicial complex X, fix an ^-norm on the sim- 
plicial chains C, [X; Z) by assigning each simplex an ar- 
bitrary nonnegative weight. (For example, every simplex 
could have weight 1, or if X is a geometric mesh in K 3 one 
could take each weight to be the length/area/volume of 
the simplex itself.) For a chain a e C n {X; Z) the Optimal 
Homologous Chain Problem (OHCP) is to find a chain c e 
C n {X; Z) homologous to a with ||c||i minimal. As there are 
many choices for c, this might seem like a hard problem. 
Indeed, consider the decision problem variant, OHCP-D: 
given X, a, and L e N is there a c homologous to a with 
||c||i<L?Weshow: 

1 .4 Theorem. The OHCP-D with integer coefficients is NP- 
complete. 

Chen and Freedman established the same result when one 
uses homology with F2 -coefficients [6]. However, with the 
addition of a simple condition on X (which often holds in 
geometric applications), Dey, Hirani, and Krishnamoorthy 
[7] have used linear programming to solve the OHCP for Z 
coefficients in polynomial time, and this will be a key tool 
here. 

We also need to consider the related Optimal Bound- 
ing Chain Problem (OBCP): given b e C„_i(X;Z), find the 
ce C n {X;Z) where d n c = foand ||c||i is minimal. (Of course, 
this is only interesting when [b] isOinH„_i(X; Z) as other- 
wise no such c exists.) If H n (X) = then the OBCP is equiv- 
alent to a related instance of the OHCP (Theorem 2.4). 
Thus one can often use the method of [7] to solve such 
problems quickly (Cor. 2.5). Conversely, we describe how 
a key construction in [1] shows that the decision prob- 
lem variant of OBCP is NP-complete in general (Theo- 
rem 4.2), and we then modify the construction to prove 
Theorem 1.4. 

To prove Theorem 1.3, we reduce the Least Area Sur- 
face problem to the OBCP via a standard desingulariza- 
tion method from 3-dimensional topology that turns an 




Figure 2.1. The surface shown is the solution to the OBCP 
for the 1 -cycle b shown in dark blue. Here X is a cube tri- 
angulated with 19,201 tetrahedra, 39,758 triangles, 24,256 
edges and 3,700 vertices, and the 1-cycle b is a Hamilto- 
nian cycle of the cube corners. The surface approximates 
Scherk's minimal surface. 

arbitrary 2-cycle into an embedded surface that is homol- 
ogous to it. This surface can be built constructively, and we 
outline in Section 6 how this gives an approach to Conjec- 
ture 1.2. As a preview to this material, we give an simple 
application of the OBCP in Section 3 to the toy problem of 
finding the shortest path between opposite sides of a tri- 
angulated square. 

2. Optimal chain problems 

For the rest of this paper, all homology will be over Z, and 
so we drop the coefficients from the notation. As in the in- 
troduction, we consider a finite simplicial complex X with 
an ( l -norm on C» (X), and recall the Optimal Homologous 
Chain Problem (OHCP): given a e C n {X), minimize ||c||i 
over all c — a + d n+ \X with x e C n+ \{X). (Here a need not 
be a cycle.) The framework of [7] is that minimizing ||c||i 
can be reformulated as minimizing some linear functional 
over the lattice points in a convex region defined by linear 
inequalities, i.e. an integer linear programming problem. 
While integer linear programming is NP-complete, when 
the matrix for the boundary map d n+ \ : C n+ \{X) — ► C n {X) 
is totally unimodular (meaning every subdeterminant is 
in {-1,0, 1}), the OHCP reduces to an ordinary linear pro- 
gramming (LP) problem, and those can be solved in poly- 
nomial time. This LP problem is the integer program with 
the integrality constraints dropped, i.e., it is the LP relax- 
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Figure 2.3. At top is the solution (green) to the OBCP for 
the 1 -cycle that is given by the knot 52 (red). The knot is 
a subcomplex of the 1- skeleton of a tetrahedral mesh of a 
cuboid shown in the lower picture. 

ation of the integer linear program. Total unimodularity 
implies that the constraint polyhedron is integral [28], and 
thus so is the solution to the linear program. 

There is a simple criterion for when d n+ \ is totally uni- 
modular. Recall that a pure subcomplex of X of dimen- 
sion k is a union of fc-simplices of X including all their 
subsimplices. We say that X is relatively torsion-free in di- 
mension n if H n {L,Lo) is torsion-free for all pure subcom- 
plexes Lo c L of dimensions n and n+ 1 respectively. Exam- 
ples include any orientable manifold of dimension n+ 1, or 
any simplicial complex that embeds in IR" +1 . It turns out 
that d„+i is totally unimodular if and only if X is relatively 
torsion-free in dimension n and so: 

2.2 Theorem ([7]). IfX is relatively torsion-free in dimen- 
sion n then the OHCPfor a e C n (X) can be solved in poly- 
nomial time. 

Turning now to the Optimal Bounding Chain Problem 
(OBCP), assume that instead we are given a lower dimen- 
sional chain b e C n -\{X) and we seek the minimum norm 
c e C n {X) whose boundary is this b. (Of course, if [b] e 
H n -\ {X) is nonzero this question is moot.) Some examples 



are shown in Figures 2.1 and 2.3. For certain X we can re- 
late these two problems: 

2.4 Theorem. Suppose a e C n [X) is such that b — d n a. 
If H n {X) = then the OHCP (for that a) is identical to 
the OBCP (for that b). Hence they have the same optimal 
solutions. 

Proof. In both problems, we seek ace C n {X) of minimal 
/ 1 -norm, so the claim is that the constraints on c are ac- 
tually the same in either case. Since H n {X) = 0, having 
b — d n c is equivalent to c = a + d n+ \x as if we have the for- 
mer then d n (c — a) = b — b = and thus there is an x with 
d n+ \x-c-a. □ 

2.5 Corollary. Suppose H n {X) = and X is relatively 
torsion-free in dimension n. Then the OBCP for b e C n -\ {X) 
can be solved in polynomial time. 

Proof. Set up the OBCP for b as an integer LP, and then 
quickly solve its LP relaxation. If b is not a boundary, then 
the feasible set is empty and we're done. Otherwise, we 
claim the solution is actually integral. By Theorem 2.4 this 
LP is identical to one for the OHCP for some (uncom- 
puted) aeC n {X) . From our discussion of the proof of The- 
orem 2.2, we know the latter LP has an integral solution as 
needed. □ 

2.6 Remark. Consider a Mobius strip embedded in IR 3 . En- 
close it in a cube and triangulate the cube with tetrahe- 
dra such that the strip is part of the 2-skeleton of the com- 
plex. Let X be the cube triangulation and b the 1-cycle car- 
ried by the topological boundary of the Mobius strip. Here 
$3 is totally unimodular, which guarantees that X is rela- 
tively torsion-free in dimension 2, as required by Cor. 2.5. 
When we solve the OBCP problem as a relaxed LP, the con- 
straints are 32 c = b. The constraint matrix 82 is not totally 
unimodular [7]. Nevertheless, Cor. 2.5 guarantees that the 
minimizer c will be integral. This does not contradict the 
result about equivalence of total unimodularity and inte- 
grality. That result says that total unimodularity is equiva- 
lent to the integrality of the constraint polyhedron for all 
right-hand sides in the polyhedral constraint equations. 
Whereas in Cor. 2.5 we are saying that the solution is in- 
tegral if the right-hand side b is in the image of 02 which 
certainly doesn't include all 1 -chains in X. 

2.7 Previous work on OBCP 

The OBCP in the trivial homology case as above has ap- 
peared in Sullivan's thesis [27] , and in the work of Grady 
[12] and that of Gortler and his coworkers [16]. When in 
addition X is an {n+ 1) -manifold, Sullivan gave a poly- 
nomial time algorithm for the OBCP based on network 
flow. This idea also appears in [16] and related work. The 
basic idea in [16] is that the cycle b of dimension n— 1 
is on the boundary of the domain with trivial H n {X). 



3 



2011/3/24 



One introduces a source and a sink and connects them 
to the centers of the top dimensional simplices on the 
boundary. These edges are given infinite capacity. The dual 
graph of the codimension-1 skeleton then forms the rest 
of the edges in the network and these edges have capac- 
ities equal to the volumes of the primal codimension-1 
faces. Then by the maxflow-mincut theorem one obtains 
a maxflow and hence an optimal chain. It is an interest- 
ing question whether such network flow methods can also 
be used to prove the more general Cor. 2.5 which only re- 
quires that X be relatively torsion-free in dimension n. 

3. A relative version of the OBCP and a toy 
problem 

In geometric applications, one often cares not about the 
specifics of the cycle b e C n -i(X) but only its homology 
class in some subcomplex Ac X. Before stating the prob- 
lem in this context, we recall the basics of relative homol- 
ogy. The relative chain groups are C n [X, A) = C n {X)/C„{A) 
which we also identify with the submodule C n (X \ A) of 
C n {X) . The boundary maps for C* {X, A) are induced from 
those of C» (X) . When C„ {X, A) is viewed as C„ (X\ A), a rel- 
ative cycle c is simply one where the support of d n c is con- 
tained in A. Thus a relative cycle gives rise to an element 
[d„c] in H n -i{A) since 3„_i o d n = 0. (This map from rel- 
ative cycles to H n -\{A)) is just the connecting homomor- 
phism in the long exact sequence of the pair [20].) We can 
now pose: 

3.1 Relative OBCP. Let A be a subcomplex ofX, and p" e 
H n -\{A). Find a relative cycle c e C„{X,A) so that[d n c] = p~ 
and\\c\\i is minimal. 

3.2 Theorem. Suppose H n {X) = and X is relatively 
torsion-free in dimension n. Then the relative OBCP for 
(5 e H n -i {A) can be solved in polynomial time. 

Proof. Suppose f> e H n -i{A) is specified by a cycle b e 
C n -\{A). First, we can quickly find a chain c e C n {X) with 
d n c-b (or determine that none exists) either by solving a 
linear system over Z [8] or by solving this instance of the 
ordinary OBCP via Corollary 2.5. 

Suppose d e C n {X) is any other chain with d n d = b' 
where b' e C n -\{A) also represents f>. Since H n {X) = 0, 
it follows that [c] = [d] in the relative homology group 
H n {X,A). Thus, solving the relative OBCP for (3 is the 
same as solving the OHCP for C n (X,A). By the discus- 
sion in Section 2, it is enough to show the relative bound- 
ary map d n+ \ : C n+ i(X, A) — - C n (X, A) is totally unimodu- 
lar. This is the case since its matrix is obtained by delet- 
ing certain rows and columns from that of the original 
3 n+l : C„ +1 (X)-C„(X). ' □ 

3.3 Toy problem 

We next give a quick application to a simple problem 
that nonetheless has all the features of our proof of Theo- 



rem 1.3 in Section 5. Let X be a simplicial complex home- 
omorphic to a square, and let L and R be subcomplexes 
corresponding to a pair of opposite sides. We will show 
that Theorem 3.2 allows us to solve the following problem 
in polynomial time: 

3.4 Problem. Find the shortest embedded simplicial path 
in the 1 -skeleton ofX joining L to R. 

See Figure 3.5 for an example. Of course, this problem can 
be solved very efficiently by a variety of algorithms. We 
give it primarily to introduce the idea of desingularization 
which is needed for Theorem 1.3. Let A - L u R. Then 
H (A) = Z® Z is generated by [v L ] and [vr] for any vertices 
Vl and vr in L and R respectively. Now f> = [vr] - [vi\ is 
in Ho (X) and so consider relative cycles c £ Q {X, A) with 
[di c] = (5. Any embedded simplicial path from L to R gives 
such a c, but of course not every such c comes from a path 
(e.g. the coefficient on some edge could be greater than 
1). However, we will show that if ||c||i is minimal then it 
does come from an embedded path. This is necessarily the 
path of minimal length, and hence we will have reduced 
Problem 3.4 to the relative OBCP for f3. Since H^X) = 0, 
Theorem 3.2 applies to let us quickly find c, and hence 
answer Problem 3.4. 

Suppose c is any relative cycle in C\ {X, A). The follow- 
ing desingularization procedure gives a collection of ori- 
ented embedded loops and arcs which give the same class 
in H\ {X, A). Along an edge e of X where c has coefficient 
w, we put | w\ strands parallel to e oriented appropriately. 
Then near each vertex v of X \ A we connect up the strands 
without introducing any crossings as shown in Figure 3.6. 
Because d\ c = at v, we can do this respecting the orienta- 
tions of the strands. That is, at the vertices, there is an out- 
going strand for every incoming one. Thus we can build a 
set of loops and arcs homologous to c. 

Now suppose c minimizes ||c||i among the relative cy- 
cles with [d\ c] = p. In the desingularization of c there must 
be at least one arc from L to R, as for instance arcs joining 
L to itself give in Hq {A) . Slightly moving this arc, we push 
it back onto the 1-skeleton X 1 of X to give another relative 
cycle d e C\ (X, A) with \d\ d] = p. The coefficient of d on 
any edge is at most that of c, and so by minimality of || c|| i 
we must have c = d . Thus c corresponds to an oriented 
path in X 1 . Moreover, this path must visit any given vertex 
at most once, since otherwise a segment of the path forms 
a closed loop which could be eliminated to reduce ||c||i. 
Thus c gives an embedded path from L to R, as claimed. 

4. NP- completeness of the OBCP and the 
OHCP 

In this section, we explain how the work of Agol, Hass, 
and Thurston [1] shows that the decision problem variant 
of OBCP is NP-complete, and then use this to prove that 
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Figure 3.5. A path which uses the fewest possible number of edges to join the two vertical sides. 



Figure 3.6. At left is part of a 1 -cycle c on a portion of a triangulated surface. At right is a desingularization of it. The result- 
ing path need not be unique, but such a path will always exist. See Section 5.3 for the higher dimensional desingularization. 



the OHCP-D is also NP-complete. Precisely, consider the 
following decision problem: 

4.1 OBCP-D. Given a simplicial complex X, a chain b e 
C„_i(X), and an L £ N, is there a chain c e C n {X) with 
d n c - b and \\c\\\ < L? 

Here the complexity is in terms of the number of simplices 
in X plus the logs of || b\\ \ and L. We will show: 

4.2 Theorem. The OBCP-D is NP-complete. 

The proof of this is essentially contained in [1], and in- 
deed they use some clever tricks to reduce more geomet- 
ric problems like Knot Genus to the more combinatorial 
OBCP-D, though they do not use the latter language ex- 
plicitly. Despite this, we include a complete proof of The- 
orem 4.2 as we need to modify the construction to prove 
Theorem 1.4. As a bonus, the simpler context of the OBCP- 
D makes the idea of [1] easier to digest for those not famil- 
iar with 3-manifold theory. 

One of two key ideas in [1] is the following construc- 
tion, which relates the OBCP-D to l-in-3 SAT, which is 
NP-complete [11, 26]. Recall that in l-in-3 SAT, we are 
given boolean variables U = {u\,...,u n } and clauses C = 
{c\, . . . , c m ], where each clause contains three literals (w; or 
its negation u{] joined by v. The question is whether there 
is a truth assignment for U so that each clause has exactly 



one true literal. We now build a 2-complex X associated 
to an instance of l-in-3 SAT by gluing together several pla- 
nar surfaces, that is, 2-spheres with (open) discs removed. 
Throughout this discussion, consult Figure 4.3 for an ex- 
ample. 

The base surface Fq has n + m + 1 boundary compo- 
nents, one labeled by the symbol K and the others by the 
elements of U U C. There is also a surface P u for each vari- 
able u e U, which has one boundary component labeled 
by u, and the others labeled by c e C for each time u (but 
not u) occurs in c. For the negation of each variable u 
there is a surface F^ with one boundary component la- 
beled u and the others labeled by the appearances of u in 
the clauses. When only u or u appears in the clauses, then 
one of the surfaces is simply a disc. 

We triangulate each surface F so that every boundary 
component consists of three 1-simplices and the number 
of 2-simplices in F is 5|dF| - 4, where \dF\ denotes the 
number of connected components of dF. For each F, we 
fix consistent orientations of its 2-simplices to create a 
relative cycle which generates H2{F,dF) = Z; we denote 
the corresponding chain in C2IF) by [F]. Now we build 
X by gluing together all boundary components with the 
same labels, in such a way that the gluings between Fq 
(where every label appears) and any F u or F Tl is orientation 
reversing; in particular, 32([Fo] + [F u \] is along the circle 
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Figure 4.3. The complex X associated to [u 2 V112V u\) A 
[u 2 v u 3 v w 3 ). 

labeled u. We let b be the 1 -cycle that corresponds to the 
boundary component of Fq that is labeled K, oriented so 
that it appears in d 2 [Fq] ■ The key lemma is: 

4.4 Lemma. The 1 -in-3 SAT instance {U, C) has a solution 
if and only if there exists deC 2 {X) with b-d 2 d and || d \\ \ < 
I + 677+ 10m. 

Proof. To start, observe that if we denote the set of literals 
by V — U U {u I u e U} then any solution d to b - d 2 d 
necessarily has the form 

d = [Fo] + ^ [-^wl f° r some kv in Z 

and the boundary of any such chain is supported on the la- 
beled circles. Since each 2-simplex lies in exactly one sur- 
face, the chains in the sum above have disjoint supports, 
and an easy calculation gives 

|| = l + 5n + 5m + \\k v \\i + 5 ^ \ k v \m v 

V€V 

where ky is the vector [k v ) and m v is the number of times 
v appears in the clauses C. Using x c ,y c ,z c to denote the 
literals that appear in a clause c, we can rewrite this as: 

||£i||i = l+5n+5fn+||*vlli+5£(|it Jtl .| + |fc yc | + |i^|) (4.5) 

ceC 

To start the proof proper, first suppose we have a so- 
lution to the l-in-3 SAT instance, and let d be the chain 
where k v is 1 or depending on whether v is true or false. 
For each variable u, one of k u and k^ is 1 and the other 
0; hence taking into account the contribution from Fq, we 



see that dd is along the circle labeled u. Now for a circle 
labeled by a clause c, as exactly one literal in c is true we 
again see that dd is along this circle. Hence dd = b, and 
(4.5) gives that ||d||i = 1 + 6n + 10m. 

Conversely, let d be a chain with dd - b and || d\\ 1 < 1 + 
6n + 10m. As dd is along the circle labeled by u, it follows 
that at least one of k u and k^ is nonzero, and hence || kv\\ 5 
n. Similarly, for each clause c at least one of k Xc ,k yc ,k Zc 
must be nonzero to ensure d has no boundary along the 
circle labeled c. Hence from (4.5) and the bound on ||d||i 
it follows that || ky II = n and each summand in the right- 
hand sum is 1 . Thus exactly one of k u and k^ is 1 and the 
other 0, and each clause has exactly one F v surface coming 
into it having nonzero weight. Therefore d corresponds to 
the needed solution to the l-in-3 SAT instance. □ 

It is now easy to prove Theorem 4.2 and then adapt this 
construction to show Theorem 1.4. 

Proof of Theorem 4.2. The OBCP-D is in NP as we can use 
the cycle c itself as the certificate. One just has to check 
that dc = b, and this matrix vector multiplication is poly- 
nomial in the number of simplices in X . Conversely, the 
OBCP-D is NP-hard since given an instance of l-in-3 SAT, 
by Lemma 4.4 there is an associated 2-complex X (made 
from 0(/2+ m) simplices) and a 1- cycle b so that a solution 
to the l-in-3 SAT problem is equivalent to finding c with 
dc- band \\c\\\ < I + 672+ 10m. □ 

Proof of Theorem 1.4. Let CX be the cone on X, which 
is a 3-dimensional simplicial complex. Let Y be the 3- 
complex obtained from CX by attaching a 2-simplex er to 
the boundary component of Fq c X c CX labeled K. For 
convenience, we use an ^-norm on C 2 (Y) so that each 2- 
simplex in X' - X u a has weight 1, but the rest each have 
weight 10 + 6n + 10m. As CX is contractible, the space Y 
is homotopy equivalent to S 2 , and hence H 2 {Y) = Z. Let 
a£C 2 {Y) generate H 2 ( Y) and have weight + 1 on a, where 
a is oriented compatibly with Fo- We claim that our in- 
stance of l-in-3 SAT has a solution if and only if a is ho- 
mologous to a cycle of weight at most 2 + 677 + 10m. Such 
a cycle would have to be confined to X' , and thus have the 
form a + d where d\s as in the proof of Lemma 4.4, proving 
our claim and hence the theorem. □ 

5. Least area surfaces bounded by a knot 

Recall from Section 1 that a basic question about a smooth 
knot K in a closed Riemannian 3-manifold Y is the mini- 
mal area of an embedded surface S with boundary K. Agol, 
Hass, and Thurston [1] considered the following discrete 
version. Take K to be a subcomplex of the 1 -skeleton of a 
triangulation of Y, where each simplex has a fixed geomet- 
ric shape corresponding to a simplex in IR 3 with rational 
edge lengths. They showed that the question of whether 
K bounds a surface of area < Aq is NP-hard. Because they 
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put no restriction on the surface involved, it is not clear 
whether this question is in NP. 

Here, we consider another discretization which is a lit- 
tle more combinatorial and will thus turn out to be NP- 
complete. For ease of exposition, let us fix that K is null- 
homologous in Y as otherwise there are no such S. We 
switch focus to the exterior of K, that is, the complement in 
Y of a small open tubular neighborhood of K. This exterior 
is a compact 3-manifold whose boundary is a torus. Let M 
be a simplicial complex triangulating the exterior, where 
each 2 -simplex has an "area" that is an arbitrary natural 
number. An orientable surface S' in Y with boundary K 
gives a properly embedded surface S - S' r\M in M whose 
boundary generates the kernel <A> of Hi [dM) — Hi (M). (A 
properly embedded surface S in M is one such that dS is in 
dM. A generator A for the kernel of Hi {dM) — Hi (M) is a 
longitude curve on dM.) Conversely, any properly embed- 
ded orientable surface S in M where [dS] = A (in Hi {dM)) 
gives a surface bounding K, after possibly adding some an- 
nuli and discs to boundary components of S to reduce the 
number of boundary components to one. 

To keep things combinatorial, we consider surfaces 
which lie in, or at least near, the 2 -skeleton of M. Initially, 
we drop the condition that the surfaces be embedded 
and consider the set & of simplicial maps /: (S,dS) — ► 
(M, dM) where S is an orientable surface with boundary 
and f*[ldS]) generates the kernel of H x {dM) — H X (M). 
(Here /„ is the map at the level of homology induced by / 
and f#, which we will use later, is the induced map at the 
level of chains.) The areas of the 2-simplices of M can now 
be used to define the area of this surface, which we denote 
Area(/). We now consider: 

5. 1 Least Spanning Area. Given Aq e N and the exterior M 
of a null-homologous knot K c Y, is there an f e 3F with 
Area(/)< A ? 

With respect to the complexity of the number of simplices 
in M and log Aq, we show 

5.2 Theorem. The Least Spanning Area Problem is NP- 
complete. 

The proof that Least Spanning Area is NP-hard is essen- 
tially the same as in [1] , and that it is in NP will follow from 
the desingularization procedure discussed below. 

5.3 Desingularization 

As in Section 3.3, a key tool is the following well-known 
procedure for turning a relative cycle ceCz [M, dM) into a 
properly embedded surface S representing the same class 
in H 2 (M, dM) . Let B be the union of small balls about each 
vertex of M\dM, and T be the union of B with even smaller 
tubes about each edge of M \ dM. For each 2-simplex a in 
M, we take oriented parallel copies of the hexagon a \ T 
according to the weight of c on a. (If some of the edges 
of a lie in dM, then a \ T may have fewer than six sides.) 



Now in the tube of T about an interior edge e of M, we 
join the adjacent hexagons to form a properly embedded 
oriented surface S in M\ B; the picture here is analogous 
to the product of Figure 3.6 with the interval, and we can 
always do this because dc = along e. The surface S meets 
the boundary of each ball B in B in a collection of simple 
closed curves in the sphere dBo. We can take a disjoint 
collection of disks in B with the same boundary as S n dB 
and add them to S. The result is a properly embedded 
surface S that is homologous to c. The way we built it, 
the surface S has the following natural decomposition as 
a simplicial complex so that the map that pushes it back 
onto the 2-skeleton of M is simplicial. In particular, we give 
S the simplicial structure where there is one triangle for 
each hexagon, one edge for each gluing of hexagons across 
the tubes of T, and one vertex for each disk added inside 
B. The desired simplicial map S ->• M just maps things to 
the corresponding simplices of M, e.g. a triangle r coming 
from a hexagon h goes to the a that h was build from. We 
summarize our discussion as: 

5.4 Lemma. Let c be a relative cycle in C2(M,dM). Then 
there is a simplicial surface S and a proper embedding S — - 
M that is arbitrarily close to a simplicial map f: S — M. 
Moreover, \\c\\i = Area(/). 

We now use this to connect the Least Spanning Area Prob- 
lem to the decision problem variant of the relative OBCP: 
given /3 e H n -\(A) and Aq e N is there a relative cycle 
c e C„(X, A) with [d„c] = jS and ||c||i < A ? 

5.5 Theorem. The Least Spanning Area Problem is equiva- 
lent to the relative OBCP-D for M and he Hi (dM). 

Combining Theorem 5.5 with Theorem 3.2 immediately 
proves Theorem 1.3, since orientable 3 -manifolds are rel- 
atively torsion-free in dimension 2. 

Proof of Theorem 5.5. Suppose that c e Cz{M,dM) solves 
the relative OBCP-D problem, i.e. [dc] — X and ||c||i < Aq. 
Then by Lemma 5.4 there is a corresponding surface / e SP 
withArea(/) = ||c||i. 

Conversely, suppose / e & with Area(/) < Aq. Con- 
sider c = f#{[S]), which is a relative cycle in C2(M,dM) and 
moreover [dc] = /* ([OS]) = A. Moreover || c\\ i < Area(/) and 
so c solves this relative OBCP-D instance. □ 

Proof of Theorem 5.2. First, the relative OBCP-D is in NP as 
we can just use c as the certificate; the sizes of the coeffi- 
cients are uniformly bounded by Aq so this is small and 
checking that [dc] = A is a polynomial time computation 
in linear algebra over Z, polynomial in the size of M. The- 
orem 5.5 now gives that Least Spanning Area is in NP. 

The argument that Least Spanning Area is NP-hard is 
essentially the same as for the original discretization of 
the Riemannian least area problem studied in [1]. Their 
proof uses a suitable 3-manifold built from the 2-complex 
of Section 4. The only modification needed here is that 
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we've set things up to require that K is null-homologous 
in Y. This can be arranged by adding a disk with weight 
10 + 6rc+ 10m to the boundary component of F labeled K, 
and also adding there a small annulus where the unglued 
boundary becomes the new K. □ 



5.6 Remark. The work of Sullivan in [27] is the closest an- 
tecedent to our Theorem 1.3. There, his motivation is rig- 
orously approximating the area of a smooth least area sur- 
face in R 3 bounding K, and the bulk [27] is showing that 
certain types of meshes necessarily contain a solution c to 
the OBCP for K whose area is within the given tolerance of 
the minimal area. While he gives a fast algorithm to find c, 
he does not insist that c gives an embedded surface. After 
all, the regularity theorems for least area surfaces guaran- 
tee that there is always such a surface of smaller area than 
c. However, once one completely discretizes the problem 
as we have done here, a priori there could be a differ- 
ence between the OBCP and the more geometric ques- 
tion about embedded surfaces. From the point of view of 
our desired application to Knot Genus as discussed in Sec- 
tion 6, it is important to have a concrete least area surface 
rather than just a homology class. 

5.7 Remark. When the ambient manifold Y is U 3 , S 3 , or 
H 3 , an important alternate approach to finding least area 
surfaces was introduced by Pinkall and Polthier in [24] , see 
also [15, 25]. They consider simplicial surfaces where the 
vertices are allowed to have arbitrary positions in Y. A dis- 
crete minimal surface is then one for which small pertur- 
bations of the vertices do not decrease total area. The pa- 
per [24] gives an algorithm that takes an initial surface So 
bounded by K and flows it toward a discrete minimal sur- 
face S. 

There are two problems with using the approach of [24] 
to solve Least Spanning Area in the restricted case ofY- 
M 3 . The first is that singularities can develop during this 
flow [24], and it seems unknown whether one still always 
ends up with a discrete minimal surface [25, §5.3] . A more 
fundamental problem is that there can be many discrete 
minimal surfaces spanning K which are not least area, and 
any flow method can get stuck on such a surface depend- 
ing on the choice of So- An extreme case is the knots of [17] 
which have infinitely many incompressible spanning sur- 
faces where no essential simple closed curve on the sur- 
face bounds a disk in M 3 \ K. Each of these surfaces should 
be isotopic to a discrete minimal surface, leading to in- 
finitely many distinct minimal spanning surfaces. 

To quickly approximate the smooth least-area surface 
spanning K in U 3 , a promising strategy is to first use The- 
orem 1.3 with respect to some mesh containing K to pro- 
duce So and then apply [24] which is not constrained by 
the initial choice of mesh. 



6. Future work 

The Knot Genus and Least Spanning Area problems have 
a very similar flavor, and hence Theorem 1.3 is compelling 
evidence for the tractability of Conjecture 1.2. However, 
these two problems are not always solved by the same 
surface — one can always cook up triangulations so that 
the least area surface is not the one of minimal genus (a 
very striking example of this is [14]). Still, for the triangu- 
lations that one encounters in practice, they should fre- 
quently be the same. Thus a natural place to look for the 
minimal genus surface is the one constructed in proving 
Theorem 1.3. 

Here, it important to emphasize how much faster the 
method of [7] is in practice compared to traditional nor- 
mal surface algorithms. Using normal surfaces, a triangu- 
lation with 30 or 35 tetrahedra is near the limit of feasi- 
ble computation, whereas [7] can handle examples with 
more than 20,000 tetrahedra (see Figure 2.1 below). Thus 
we should be able to work with reasonably fine triangu- 
lations of M, which could increase the chance that the 
least area surface is minimal genus. For instance, from the 
Thurston/Perelman point of view, one could take some 
combinatorial approximation to a hyperbolic metric on 
M. Unfortunately, doing so will not completely eliminate 
the issue of least area surfaces not having minimal genus; 
in the smooth category, a folk theorem gives a hyperbolic 
3-manifold with a homology class whose minimal area 
representative is compressible (cf. [13]). 

However, going to a large number of tetrahedra raises a 
different issue: we need to know the genus of the surface S 
constructed from the ( l -minimal cycle c. In fact, the sur- 
face S is not unique, as there can be many choices for how 
to pair up the triangles near the edges of M (cf. Figure 3.6), 
sometimes resulting in exponentially many such surfaces 
in terms of || c|| i . Moreover, these choices affect how many 
disks we add near the vertices of M to complete the sur- 
face, and thus affect the Euler characteristic and hence the 
genus. This leads to the natural question: Can one quickly 
determine min (~x(S)) over all surfaces S resulting from c? 
It might also be the case that the classes c that one finds in 
practice don't have many surfaces associated to them, and 
we will run computer experiments on this. 
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